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Abstract 



i-C I In the framework of unitary Glauber approximation for particle- 

1^ atom scattering, we have received an improvement of the Moliere-Fano 

Mk theory (M-F theory) on the basis of reconstruction of the generalized 

r~| ! optical theorem in it. We have estimated the relative unitary and the 

Coulomb corrections to the parameters of the M-F theory and studied 
^ ■ the Z-dependence of these corrections. We showed that the absolute 

^ I and relative inaccuracies of the Moliere's theory in determining the 

screening angle for heavy atoms of the target material {Z ~ 80) are 
Q ■ about 20% and 30%, respectively. 

|>; PACS: ll.80.La, 11.80.Fv, 32.80.Wr 
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• r-j 1 Introduction 

>< 

c^ I The Moliere-Fano multiple scattering theory of charged particles' [1-3] is the 

most used tool for the multiple scattering effects accounting at experimental 
data processing. The experiment DIRAC [U |5] like many others |6j (MuScat 
[7], MUCOOL [8] experiments, etc.) meets the problem of the excluding of 
multiple scattering effects in matter from obtained data. 

The standard theory of multiple scattering [H El E] , proposed by Moliere 
[H 12] and Fano [3] , and some its modifications [7-12] are used for this aim. 
The modifications, developed in Refs. [7-9], are motivated by experiments 
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[3 E] and are connected with the inclusion to Mohere theory of analogues 
of Fano corrections, as well as the determining their range of applicability 
[7~10]. In Ref. [H] is presented a modified transport equation whose solu- 
tion is applicable over the range of angles, from to 180°. In [T2] results 
of experiments [T3] are qualitatively explained within the framework of the 
theory allowing for pair correlations in the spatial distribution of scatterers. 
Estimation of the theory accuracy is of especial importance in the case of 
DIRAC experiment for it's high angular resolution. One possible source of 
the inaccuracy of the M-F theory is use in [1-3] an approximate expression 
for the amplitude of target-elastic scattering of charged particle by atom 
which violates the generalized optical theorem 

k k 

47r 47r 

or, in other words, unitarity condition. Another possible source of inaccuracy 
is using in calculations an approximate relation for the exact and the Born 
values of the screening angle (Xa) 



obtained in the original paper Moliere ^Lj- Therefore, the problem of esti- 
mating the M-F theory accuracy and improvement of this theory becomes 
important. 

In the present work, we have estimated the relative unitary corrections to 
the parameters of the Moliere-Fano theory, resulting from a reconstruction 
of the unitarity in the particle-atom scattering theory, and found that they 
are of order Za^. We have also received rigorous relations between the exact 
and first-order parameters of the Moliere multiple scattering theory, instead 
of the approximate one obtained in the original paper by Moliere. Addition- 
ally, we evaluated absolute and relative accuracies of the Moliere theory in 
determining the screening angle and have concluded that for heavy atoms of 
the target material they become significant. 

The paper is organized as follows. In Section 2, we consider the approx- 
imations of the M-F theory. In Section 3, we obtain the analytical and 
numerical results for the unitary and the Coulomb corrections to the param- 
eters of the M-F theory. In Conclusion, we briefly summarize our results. 
The some preliminary results of this work were announced in Refs. |14j . 

The given work is dedicated to the 
70th Anniversary of Alexander Tarasov 



2 Approximations of the M— F theory 

2.1 Small-angle approximation 

Let all scattering angles are small ^ ^ 1 so that sin 6^6, and the scattering 
problem is equivalent to diffusion in the plane of 9. Now let cr^iix) be the 
elastic differential cross section for the single scattering into the angular 
interval x = — 9', and WM{0,t)9d9 is the number of scattered particles 
in the interval d9 after traversing a target thickness t. Then the transport 
equation is 

^^gf'^^ = -noW,,{9,t) I a,^{x)d'x + no I WUO-X,tKiix)d'x, (1) 

where no = cNa/A is the number of the scattering atoms per cm^, Na is the 
Avogadro number, and A is the atomic weight of the target atoms. 

Moliere solved this equation for the determination of the spatial-angle 
distribution function WM{9,t) by the Fourier-Bessel method and gets to her 
a general expression 



W,,{9,t) = j M9ri)g{T],t)'ndr], (2) 



in which 

^(r/,t)=exp[iV(r/,t)-iVo(0,t)], (3) 

9 is the polar angle of the track of a scattering particle, measured with respect 
to the initial direction 2;; 77 is the Fourier transform variable corresponding 
to 9; and the Bessel function Jq of order is an approximate form for the 
Legendre polynomial appropriate to small scattering angles [21 [15] . 
In the notation of Moliere 



N{r],t) = 2Tmot / a^i{x)Mxv)xdx, (4) 



and A^o is the value of (jl]) for 77 = 0, i.e. the total number of collisions 

00 
NoiO,t)=2nnotfa^i{x)xdx- (5) 



The magnitude oi Nq — N is much smaller than A'"o for values 1], which are 
important, and may be called the 'effective number of collisions'. 
Inserting Eqs. ©-([Sj) back into ([2]), we have 







-27mot / a,;(x)xc?x[l - Mxv)] 





(6) 



This equation is exact for any scattering law, provided only the angles are 
small compared with a radian. 

At g{ri, 0) = 1 for all rj the expressions ([2])-([S]) can be rewritten as follows: 

oo 

W,,{e,t) = y Jo(er/)e-"°*'3-('')r7rfr/, (7) 



where 

CX) 

QM = 2nJaM)[l - Mxv)]xdx- (8) 



This result is mathematically identical with result of Snyder and Scott for 
the distribution of projected angles [T6] . 

2.2 Approximate solution of the transport equation 

One of the most important results of the Moliere's theory is that the scat- 
tering is described by a single parameter, the so-called 'screening angle' (xa 

or Xa)- 



X'a = VTierXa = [exp (C, - 0.5)] Xa ~ 1.080 Xa, (9) 

where Ce = 0.57721 is the Euler's constant. 

More precisely, the angular distribution WM{0)6d6 depends only on the 
ratio of the 'characteristic angle' Xci which describes the foil thickness, to the 
'screening angle' which describes the scattering atom: 

h^^=\n(^\ =\^{—\ +l-2C^~lniVo. (10) 

The screening angle Xa can be determined approximately by the relation 

Xl ^ Xo(l-13 + 3.76a2) = {Xa? (l + 3.340^) (11) 



with the so-called 'Born parameter' 

_ Ze^ _ Za 



(12) 



The second term in flTTl) represents the deviation from the Born approxi- 
mation. If the value of this term is 0, the value of the screening angle is 

The angle Xo is defined by 



Xo = 1.13' 



^1/3^ ^1/3 



niea 



137p 0.885p ' 



(13) 



where p = meV is momentum of incident particle and v is its velocity in the 
laboratory frame. 

The characteristic angle is defined as 



Xl = 4.n,i (1^) 



(14) 



Its physical meaning is that the total probability of single scattering through 
an angle greater than Xc is exactly one. 

Putting XcV = y cind setting 9/xc = u, for the most important values of 
7] of order l/xc we get Moliere's transformed equation 



WM{0)6d6 = udu I ydyJo{uy)exp <^ — — 



&e.-ln(^ 



(15) 



which is very much simpler in form than (|6]). 

In order to obtain a result valid for large all angles, Moliere defines the 
new parameter B by the transcendental equation 



B-\nB = b^i. 



The angular distribution function can be written then as 

OD 

WU0,B) = l^JydyMey)e-y'/'exp 



AB V 4 



(16) 



(17) 



The presented expansion method is to consider the term [?/^ln(?/^/4)]/4i? as 
a small parameter. Then the Wm is expanded in a power series in 1/ B\ 



1 1 



^-(^'^) = E;;y^^"(^'^) 



n=0 
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with 



oo 

2 /„,2\ -[■"■ 



Wn{0,t) = i- fydyjJiy]e-y'" 



02 J - '-^ -yQ 





y ^^(y 



4 V 4 



(19) 



P = xlB = A^not(—\ B{t). 
\pv J 



This method is vahd for B > 4.5 and 6'^ < 1. 

The first function Wo{6,t) have simple analytical form: 

Woie,t) = lexp(-ty (20) 



^2 ~ tint. (21) 

i— > oo 



For small angles, i.e. 9/6 = 6/(xcy/B) = O less than about 2, the Gaussian 
fl20|) is the dominant term. In this region, Wi{6,t) is in general less than 
Wo{6,t), so that the corrections to the Gaussian is of order of 1/B, i.e. of 
order of 10%. 

2.3 Born approximation 

On the one hand, Moliere writes the elastic Born cross section for the fast 
charged particle scattering in the atomic field as follows: 

<{X) = cr^x) [l - ^^^ =<{x)q!i{x). (22) 

For angles x small compared with a radian the exact Rutherford formula has 
a simple approximation: 

"" «5(X). (24) 
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vmot X 

Here, F^ is the atomic form factor and q^iix) is the ratio of actual to the 
Rutherford scattering cross sections in the Born approximation. 



Then the screening angle xf in the Born approximation one can represent 
via F4 or q^i{x) by the equations 



In (Xf ) 



hm 



hm 

i—^oo 



Lo 






^IM^ + l_lnc 
X 2 



with an angle <;" such as 



Xo < ^ < 1/^ ~ Xc 



(25) 
(26) 



(27) 



and the angle Xo ~ rUeaZ^/^ /p. 

Moliere approximation for the Thomas — Fermi form factor Frp^p{q) with 
momentum transfer q can be written as 



Fr-Aqr = E 



cAl 



-q' + X. 



2 ' 



(2^ 



in which 



ci = 0.35, C2 = 0.55, C3 = 0.10, 

Ai = 0.30A, A2 = 4Ai, A3 = 5A2. 

In the case where the Born parameter a = 0, the equation (I25p for the 
screening angle can be evaluated directly, using the facts that g(0) = and 
lim q{() = 1. Then with use of fl22l) and fl28l) can also be obtained the 

following approximation for (^Xa) ID [IS] : 



A 



iXa) = [exp(C, - 0.5)] -A = VUU xo^, 



(29) 



where A = nieaZ^^^ /0.885. We should note that in Refs. [U dH] admitted a 
misprint, namely, the factor A = 1.0825 in Eq. ( |29|) should be replaced by 
A = 1.065 = Vtl3. 

On the other hand, Moliere writes the nonrelativistic Born cross section 
in the form 



c^!iix) = k' 



00 
y"prfpJo(2A;psin|)K(p) 



(30) 



where the Born phase shift is given in units oi h = c = 1 hj 

OO CO 



-OO 



Here, k is the wave number of the incident particle, the variable p corresponds 
to the impact parameter of the collision, and U\{r) is the screened Coulomb 
potential of the target atom 

Ux(r) = ±Z-M\r) (32) 

r 

with a screening function A(Ar). The screening radius A~^ = r^^ is frequently 
taken to be the Thomas — Fermi (T-F) radius Vgc = 0.885 /irieaZ^''^. 

Snyder and Scott [16] have used the simplest form of (l32l) with a pure 
exponential factor A(Ar) = exp(Ar): 

Ux{r) = ±Z-e''. (33) 

r 

Moliere in his detailed study of single scattering [1] approximated the T-F 
screening function by a sum of three exponentials: 

A(Ar) ~ 0.1e-^^" + 0.55e-i'2^" + 0.35e-°-3^". (34) 

For the Born target-elastic single cross section have a place the following 
relations: 



3^(0) = ^"St^-^"^. (36) 

The Born approximation result for the amplitude fei of elastic scattering of 
charged particle by atom reads 



Ue) =ik j Upq)[l - e^*«(^^]prfp, (37) 



q = ke. 



2.4 Approximate relation for the quantities Xa and x. 



In order to obtain a result valid for large a (i.e. not restricted to the Born 
approximation), Moliere uses in his calculations of the screening angle a 
WKB technique. 

The exact formulas for the WKB differential cross section o"ei(x) ^ind 
corresponding qei{x) ^^ terms of an integral are given in Moliere's paper pQ 
as follows: 



(^elix) = k"^ 



oo 

pdp Jo{kxp) < 1 - exp [i^t 



(St 



Qelix) 



4a2 



pdpJo{kxp)\ 1 - exp [i$Af(p)] 



with the phase shift given by 



m 



krir) — k 



dz. 



(40) 



where kr{r) is the relativistic wavenumber for the particle at a distance r 
from the nucleus, and the quantity p is seen to be impact parameter of the 
trajectory or 'ray'. As before, k is the initial or asymptotic value of the 
wavenumber. 

When kr{r) is expanded as a series of powers of Ux{r)/k, then the first- 
degree term yields the same expression for ^m{p) as Eq. fl^ . The Born 
approximation for (1551) is obtained by expanding the exponential in ( 155]) to 
first order in parameter a (TT2|) . 

Relations (^^ and (l52l) between the quantities cr^i{x): Qeiix)^ and Xa 
remain valid also for o"e;(x), '?e/(x)) and Xa- 

Despite the fact that the formulas ( 1551) . ( 15^ are exact, evaluation of these 
quantities was carried out by Moliere only approximately. 

To estimate (l39|) . Moliere used the first-order Born shift (J3T]) with (l32|) 
and (IM|) . what is good only to terms of first order in a, and found 



Qelix) 



1 



4ia{l — ia) 

[xIxqY 
1 



0.81 + 2.21 



- ^ [i){ia)] 



la 2ia 2xo 



(41) 



Here, if) is the so-called 'digamma function', i.e. the logarithmic derivative 
of the F-function ^(x) = d\n.T{x)/dx. 

He has fitted a simple formula to the function 3f? ['?/'(«a)] in Eq. (HTj) 



3fJ[^(^a)]^ilg(^a^ + ^ + 0.13^. (42) 

Using f H2|) in Eq. ( 14T|) and expanding with neglect of higher orders in a? and 
(x/Xo)"^ he got 



8.85 



oj 



'i^oqnq 2,^ 7.2-10-%/xor 

1 + 2.303 a IgT^ — 2 /q , m q^ 

[a* + a"'/3 + 0.13) 



(43) 



Moliere has calculated Qeiix) fo^' different values of a and as a result has 
devised an interpolation scheme, based on a linear relation between {x/Xof 
and a? for fixed g^;: 

(X/Xo)' ~ A + «'5<?- (44) 

Calculating the screening angle defined by 

(leAx)dx 



- In (Xa) = - + lim 



X 



— In^ 



1 



and assuming a linear relation between x^ and a^, Moliere writes finally the 
following interpolating formula: 



Xa ~ XoVl.l3 + 3.76a2. (46) 

Critical remarks to his conclusion see in the review of V.T. Scott ITGl. 



2.5 Fano approximation 

To estimate a contribution of incoherent scattering on atomic electrons is 
often [21 [31 [151 [E] replaced the squared nuclear charge Z'^ with the sum of 
the squares of the nuclear and electronic charges Z{Z + 1) in basic relations 
(for differential cross-section, some parameters of the theory, etc.). 

This procedure would be accurate if the single-scattering cross sections 
were the same for nucleus and electron targets. Besides, the actual cross 
sections are different at small and large angles. Fano modified the multiple 
scattering theory taking into account above differences. 
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For this purpose Fano separates the elastic and inelastic contributions to 
the cross sections 



O-(X) = (^elix) +(^in{x)- 



(47) 



For the inelastic components of the single scattering differential cross 
sections, the Fano approximation reads 



daj' 



da. 



B 



dn 



dn 



(4^ 



Since the Born single-scattering amplitudes are pure real, the generalized 
optical theorem cannot be used to calculate the total cross section in the 
framework of this approximation. 

Fano sets the task of comparing the contribution of aj^ix) to the exponent 
of the Goudsmit — Saunderson distributioio [ISJ for total scattering angle: 



w{e,t) 



2-E(' + 5) 



Pi{e) exp -no t / a^(x) smxdx[l - Pi{x)] , (49) 



where Pi is the Legendre polynomial. If we replace in fH9|) the sum over / 
by an integral over r], (/ + i) by t], Pi by the well-known formula Pi{0) = 
Jo( (^ + ^) 6*), and sinx by x, the expression fll^ goes over into small-angle 
distribution ([H]) of Moliere and Lewis. 

To achieve the mentioned goal in the small-angle approximation, we de- 
termine in accordance with p^ and [9j similarly to ( 123|) and (125|) the corre- 
sponding expressions for the inelastic cross section 



<ix) = (^''{x)qUx) 



Xc 



Arrriot Z [1 — cosx)^ x'^ 



itnix) 



(50) 



Xc 



vmot Z x^ 
and the 'inelastic cut-off angle' XiL 

C 



C(x) 



- In (xfn) = lim 



X 2 



(51) 



(52) 



^ The Goudsmit — Saunderson theory is valid for any angle, small or large, and do not 
assume any special form for the differential scattering cross section. 
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Then, with the use of (I2S]) and ( IFTj) . angular distribution (IT71) can be 
rewritten as follows 



oo 

W,,.,{e,B) = i^ j ydyJo{ey)e-y"/'exp{Y^i + Yj (53) 



with 

oo 

c 

where the parameter B is defined by equation 

B-\nB = b^i + b,^, (55) 

in which 

^,^.n(|)%i-2C.. ..„.^,„(|)^ ,56) 



^.jn 



Numerical estimation of the — Mj„ = — In (Xin) value from the T-F model 
yields (— Mj„)j,_^ = 5.8 for all Z. The value of — Mjn should not vary greatly 
from one material to another. 

For sufficiently large angles, with the use of exact Rutherford formulas 
([23]) and (15U]) . the correct angular distribution W{9,t) may be estimated 
according to the formula 

W,UO,t) = Wie,t) [af_,(x)/a^(x)], (57) 

as suggested Bethe [U] and Fano |3]. 

3 Improvement of the M— F theory 

3.1 Glauber approximation 

The Glauber approximation [12] for the amplitude Fif[q) of the multiple 
scattering at small angles is given by 

F^f{q} ='^ld'p exp(^ gt) T,f{p}, (58) 
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where Tif{p) is so-called 'profile function'. 

We may formulate the problem in a general way by considering the scat- 
tering of a pointlike projectile on a system of Z constituents with the coordi- 
nates fi, r2, . . . jfz and the projections on the plane of the impact parameter 
Si, 5*2, •• • ,Sz- Then we may write the total phase shift as a sum of the form 

z 
X(p,xi, X2,... ,Xz) = ^Xkip- Sk). (59) 

fc=i 

If we introduce the configuration space of the wave functions ipi and ipf for 
the initial and final constituents states, we may write the profile function as 

z 

r./(p) = l['^'nrf{{n})^P^{{rk})r{p,{sk}) (60) 

"^ fe=i 

with the interaction operator 

r(p, {sk}) = 1 - exp[z<l>(p, {sk})] (61) 

and the phase-shift function 

z 

$(p, {4}) = Zxip) -Y^xip- 4). (62) 

fc=l 

When the interaction is due to a potential V{r), the phase function x{p) is 
given by 



oo 

-- f V Up^ + zA dz. 



X{P) = -- V{^/^+^)dz. (63) 



with the potential of an individual constituent's 

V(r) = ± lim - e-^'', X ~ meaZ^/^. (64) 

The multiple-scattering amplitude Fif{q) fl58|) is normalized by the rela- 
tions: 

ATiQFuiO) = a{i)tot, \Fif{q)\^ = da^f/dq^, (65) 

where 

(^{i)tot = <y{i)ei + cr{i)in, (Tif= hdaif/dgrl'^d'^q, (66) 

13 



a(.) 



tot 



/ 



(67) 



In terms of e**, where the phase-shift function $ = $(p, {sfc}) is given by 
62|) . the cross sections a{i)tot, o'{'i)ei, <^{i)in could be written as follows: 



a[t)tot = 2^ I (l-(e**))rfV, 



(68) 



a(z) 



el 



l-23fj: 



(e^*) + |(e^^)|')rfV 



(69) 



crli 



,j<j> 



rfV- 



(70) 



The brackets ( e'* y signify that an overage is be taken over all configurations 
of the target constituents in i-th state. 



3.2 Reconstruction of unitary conditions 

To reduce the many-body problem to the consideration of an effective one- 
body problem and to establish the relationship between the Glauber and the 
Moliere theories, we introduce the abbreviation 

(e^*) = e'* (71) 

and consider for the effective ('optical') phase shift function $(p) the expan- 



sion 



°o ^n-1 



^(/^) = E^* 



(72) 



71=1 



where 



$, 



$! = ($), $2 = (($-$1)2), 

= ((<l>-$i)3),...,$„~Za7/3. 



(73) 



The first order of the expression for $(p) is simply the average of the 
function $(p, {sk}) and correspond to the first-order Born approximation. 
The second order term of $(p) is purely absorptive (i.e., positive imaginary) 
in character and is equal in order of magnitude to the Za^ / [5. 
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When in the expansion (172]) 

'^tailip) = V ^-<^n < 1, 



°o ^„_l 



n=3 



n! 



(74) 



it seems natural to neglect the ^taii{p) and consider the following approxi- 
mation: 



$(p)^ $!(/?) + -<l>2(p), 



(75) 



in which we let $i(p) = $a/(p) and $2(p) = 2$j„(p). The last term corre- 
spond to target-inelastic scattering, also referred to as incoherent scattering 
or 'antiscreening' [5]. 

This leads to the following improvement of the Moliere-Fano theory: 



'^MiP)^'^^AP)+^'^^n{p) 



with 



2<l>,„(p) = lim Z { / \xxip- TrW Qi^d'r - 

A-i>0 



Xx{p-rT)0{r)d^r 



(76) 



,(77) 



where 



Xxip) 



II V, (v/7 + ^) dz, V, 



a 
r 



-Xr 



{n 



Qif) = ^}{r}^i{r). 
For the cross sections 

(^tot = {cr{i)tot) , CTin = {(j{i)in) , a^ = {cr{i)ei) 
takes place the unitary condition 

k k 



with 



fei{e)=Fu{q}, ^Fu{0) = ^J2 J \^^f(^\' ^^^ 



(79) 



^0) 



^2) 
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ri/(p) = l-exp(-2$,„), (84) 

Using (ISTj) . we found for the cross sections the following expressions: 



atot = 47r / (1 - cos$,,(p) e-*-^'^) prfp, 
ain = 2nf{l-e-^''^"^^)pdp, 



15) 



(86) 



aez = 271 f {1-2 cos $,,(p) e-*-^'^ + e-^*-^'^) pdp. (87) 

3.3 Unitary corrections to the Born approximation 

With the use of the evaluation formula 

[2^Up)]d^p-Zay/3 (88) 

and the exact contributions have been calculated in [20], we obtain the fol- 
lowing unitary relative correction 6unit to the first-order Born cross section 
of the inelastic scattering a^^: 

Sunu = ^ = ^^^^^ = ^ - 1 ~ Za'/i^ (89) 



in in in 



with 

al = (^j ^in{p)dPp\. (90) 

The corresponding angular distribution reads 

oo 

Win{e) = 2n f Joier])e-'^'-^'^^r] dr], (91) 
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QUv)=27r / aU0)[l-M0v)]9d9. 



Inserting Eq. fl92|) back into flOTj) . we get the equation of the form: 



Winie) = 27T / r]dr]Jo{e7])exp 



With the use of 



-2tt / a 



\i\ni jni 



de'[i-Me'r^)] 



r]J,{9r])Jo{9'r])dr] = -6{e-9') = 



and 



(92) 



(93) 



(94) 



c^^^ Joi^v) = 2a" 



r(i) 
r(o) 



0, dr]riM9r])Jo{e'ri) = 0, (95) 





according to [21], the integration of Eq. f l93|l yields the following result: 



W,n{e) = -{2nf r]Mer])Jo{e'r])d7]-a^n{e')e'de' = -{27r)^aU0)- (96) 



In Eqs. (194|) and (l95l) S is the Dirac delta function and F is the Euler Gamma 
function. 

Finally, taking into account the relations (1891) and (I96p . we can estimate 
the unitary correction to the angular distribution function (PT]) : 

*.„.. = ^^ji^ = ^^^ - 1 = ;;ji^ - 1 = ^ - 1 ~ Za //S. (97) 



a, 



3.4 Rigorous relations between exact and Born results 

To obtain an exact correction to the Born screening angle (xa) in small- 
angle approximation, we will carry out our analytical calculation in terms of 
the function Qeiijf) dHD: 



QeK^) = 27r / aeKx)[l - Mxv)]xdx = U'pl- co4A$(p,r/)] , (98 



where the phase shift is determined by the equation 
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A$(p,r/) = $(p+)-$(p_), p± = p±r//2p. 
Substituting the expression for the cross section 



(^elix) 



xl 



into fl98l). we rewrite it in the form: 



Qeiix) 



(99) 



(100) 



uotQeiiv) = 2xc / [1 - Mxv)] qei{x)x '^dx- 



(101) 



For the important values of rj of order 1/xc or less, it is possible to split the 
last integral at the angle <;" (l27j) into two integrals: 



I{r]) = / [1 - '^o(x^)] qei{x)x dx 



/ [1 - ^o(x^)] qeiix)x'^dx + /[I - ^o(x^)] qei{x)x'^dx 



hi^ri) + hi^^r]) . 



(102) 



For the part from to <;", we may write 1 — Jo{xv) = X^'7^/4, and the integral 
Ji reduces to a universal one, independently of rj: 



A('?^) = ^ / 1ei{x)dx/X- 



(103) 



For the part from ^ to infinity, the quantity qei{x) can be replaced by unity, 
and the integral I2 can be integrated by parts. This leads to the following 
result for I2: 

,2 



hi^^v'. 



4 



l-ln(gr7) + ln2-C^ + 0(^r/) 



(104) 



Integrating (I103p with the account fH51) . substituting obtained solutions back 
into (llOip . and also using definition 

ln(xc/Xa)' + l-2C, = ln(xc/Xa)', 
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as a result for Q^iir]), we obtain: 

iXcVf 



Qeiiv) 



ln(^'l-ln^^^ 



2nnt 



iXcV)\^ (vHxaf 



Xa 



(105) 



2nQt 



Finally, considering definition of Xc flT^ - we can represent Qeiiv) by the 
following expression: 



Q.,M = -2.(|^)',= l„^"^(^-"'^' 



(106) 



Then, the screening angle Xa can be determined via Qdijf) by a linear equa- 
tion: 



-ln(xa' = l"(^) + 



21X1^' 



Jp) 



Qeiiv)- 



(107) 



Let us present the quantity Qeiiv) i^ the form: 

Qeiiv) =Q!liv)-^Cc[Qeliv)]- 



(10^ 



Then, on the one hand, using (I106p . the difference AcdQeiiv)] between the 
Born approximate Qeiiv) and exact in the Born parameter values of the 
quantity Qez(^) '^^^ be reduced to a difference in the quantities In (x^) and 



AaciQeliv)]^ Qeiiv) -Qeiiv) 



(109) 



""''il^ 



ln(x:)-lnU:)*'|^4^^'||^| Aec[ln(x:)l. 



On the other hand, this difference can be reduced to a difference Agez(x) 

Qelix) ~ qelix)- 



AaciQeiiv)] = 2nlxdxAaM)[l - Mxv)] = ^ J %Aq,iix)[l - Mxv)], 
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and for (15^ with (ESD, as a result of integration, ( IllUp becomes 



r)^0 \ PP 



^H'f)4<-'fi-^'W 



47rr7' 



Za 

Jp 



ipi l + i 



Za 

T 



+ a 



where 



3fJ [V^ (1 + ia)] = 3f? [V^ (1 - m)] = 3fJ [^ (ia)] = 3fJ [^ (-m)] 



-C^. + a^ 5^ 



n=l 



n(n^ + a?) 



-C^ + /(a) 



(110) 



(111) 



-oo < a < oo, 

2 V^oo r^/^2 I ^2\1^1 



V^(l) = —Cjs, and /(a) = a Yl^=i V^i.tT^ + " )] is an 'universal function of 
a = Za/ (3 \ 

As consequence, we can get the following rigorous relations between the 
quantities In (x^) and In {xaY- 

ln(x:)-ln(x:)" = 5R[^(l + za)-^(l)], (112) 

Acc[ln(x:)] = a'YT [n{n' + a'rK (113) 

•^— »n=l 



We point out that the relations (IllOp . (11121) . and (I113p are independent 
on the form of electron distribution in atom and are valid for any atomic 
model. 

From (IllUp also follows an expression for the correction to the exponent 
of Qj. Since ln[5f(r7)] = —riotQ, we get: 



Aec[ln^(r7)] = ln[^(r^)] - ln[^-(r^)] 
= 4:Tir]^nA-^ 1 f{a 



(114) 



For the specified value of rf = l/xl, using the definition of Xc (fT4|) . we can 
evaluate this correction: 



A.c[i^,ix.)]-'-^(^^\m^m. 



(115) 
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The formulas for the so-called 'Coulomb correction' (CC), defined as a 
difference between the exact and the Born approximate results, are known 
as the Bethe — Bloch formulas for the ionization losses [22] and the formulas 
for the Bethe — Heitler cross section of bremsstrahlung [23j. 

The similar expression was found for the total cross section of the Coulomb 
interaction of compact hadronic atoms with ordinary target atoms [21] . The 
more complicate formal expression for CC was derived by I. 0verb0 in |25] . 

Were also obtained CC to the cross sections of the pair production in 
nuclear collisions [26| |27] . as well as to the solutions of the Dirac and Klein- 
Gordon equations [28]. Specificity of the expressions received in the present 
work is that they define CC not to the cross sections of the processes, but to 
screening angle and an exponential part of the distribution function. 



3.5 Relative Coulomb corrections to the Born approx- 
imation 



Let us write flllSp in the following way 

(Xa) = {XaY ^W 



At) 



:ii6) 



The exact (Coulomb) relative corrections to the Born screening angle (xa) 
can be then represented as follows 



^Cc{Xa) 



X^-(x:)'_ AU) 



ixLY 



ixLY 



exp 



:^(t)j-- 



(117) 



The relative CC to the exponent g^{r]) at 1]"^ = l/Xc? ^s follows from (IllSp . is 
also determined by this quantity: Scc{Xa) = ^cc [diXc)]- Moreover, because 



oo 

AWixc, t) = W,,-W^ = j JomAg{xc)vdv, 



:ii8i 



accounting J^ drjrj Jo{9ri) = 0, we get 



5cc [WmKXc, t)\ = jjy^_ — j^ = _^,_ , = exp 



W^iXct),, g^ixc) 



Thus, 



/ 



/Za 



Jcc 



^cc{Xa) = ^cc [g{Xc)] = Sec [WM{Xc,t)] 



(119) 
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To compare the second-order relative corrections to the first-order results, 
which correspond to Eqs. flll6p and f ITT]) . respectively, we first present these 
equations in the following approximate form: 



\Xa) '^ \Xa) 



1 + 1.204 




(120) 



yXa) ^~~^ \Xa) 



1 + 1.670 



Za 

T 



+ 




:i2ii 



The last expression follows from 




3.34 



Za 

T 



Xa 



1 + 1.670 



Za 



+ 



where Xa = (Xa)/l-080 and xf = (Xa)7l-080. Then 
becomes 




mx'a) 



{x'aY 



1.204,^ 



^Ht) 




and (IT^ 



'1221 



(123) 



Additionally, in order to assess the accuracy of the Moliere theory in 
determining the screening angle xLi ^^ ^^^ define the absolute and relative 
differences between the values of 5m \x'a) ^^nd 5cc{x'a) ^^ 'well as between the 
5m (Xa) and 5cc {x'a) by the relations 



'JCCM \Xa) 
5cCM [x'a) 



A 



(2) 



{x'a) 



^(2) 



^CCM (Xa) 



r(2) _ .(2) 
^M ^CC 

'~>M '~>CC 



'Jm 



6 



(2) 



1- 



fJcc 

'Jm 

Sec 



(124) 
(125) 



The numerical estimation of the second-order corrections shows that the 

(2) 

values of the relative corrections 6cc for heavy target atoms with Z ~ 80 does 
reach 40%. From it is also obvious that with the rise in the nuclear charge 
the absolute inaccuracy AccM{Xa) of the Moliere's theory in determining the 
screening angle xL increases to approximately 16%, and the corresponding 
relative error 6ccM{Xa) '^o^s not depend on Z and is about 28%. 
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To estimate the exact ('Coulomb') corrections according to the formulas 
(11131) and (11171) . we must first calculate the values of the function /(a) = 

^[^{1 + ta)] +Ce. 

From the digamma series |29j 

^ oo 

^(l + a) = l_C,--^ + V(-l)"[C(ri-l)]a"-\ |a| < 1, (126) 

1 + a ^-^ 



n=2 



where ( is the Riemann zeta function, leads the corresponding power series 
ior ^[iIj{1 + ia)] = ^[^p{ia) 



^ oo 

^[iPiia)] =1-C^-^— ^ + ^(-l)"+i[C(2n+l)]a2", |a| < 2, (127) 



n=l 
2 Y^oo r / 2 I ^2M-1 



and the function /(a) = a Xl^i ['^('^ + o- )] can be represented as follows: 

^ oo 

^^^) = l-T^^ + E(-l)"^'[^(2n+l)]a2", |a|<2, (128) 



+ 
1 



n=l 



l + a2 



0.2021 a^ - 0.0369 a^ + 0.0083 a^ - . . .(129) 



The calculation results for Z dependence of (11281) . the relative CC and 
their difference with the Moliere corrections, as well as the unitary corrections 
(197|) at /3 = 1 are presented in Table 1. 

Table 1. The Z dependence of the corrections defined by Eqs. 

dnzD, dnsD, dnzD, and ([nsD 



M 


Z 


Za 


10 6 unit 


fiZa) 


^cc 


^CCM 


(>CCM 


Be 


4 


0.029 


0.002 


0.001 


0.001 


0.000 


0.286 


Al 


13 


0.094 


0.007 


0.011 


0.011 


0.004 


0.280 


Ti 


22 


0.160 


0.012 


0.031 


0.031 


0.012 


0.280 


Ni 


28 


0.204 


0.015 


0.049 


0.050 


0.020 


0.287 


Mo 


42 


0.307 


0.022 


0.105 


0.110 


0.047 


0.297 


Sn 


50 


0.365 


0.027 


0.144 


0.154 


0.068 


0.306 


Ta 


73 


0.533 


0.039 


0.276 


0.318 


0.157 


0.330 


Pt 


78 


0.569 


0.041 


0.307 


0.359 


0.182 


0.336 


Au 


79 


0.577 


0.042 


0.313 


0.367 


0.187 


0.337 


Pb 


82 


0.598 


0.044 


0.332 


0.393 


0.205 


0.342 
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The Table 1 shows that while the values of relative unitary corrections 
5unit for heavy atoms of the target material does not reach 0.5%, the maxi- 
mum value of the difference A^c [in (Xa)] = fi^a) is close to 33%, and the 
magnitude of the corresponding relative Coulomb correction 6cc is reached 
for Z ~ 80 the value of order of 40%. 

From it is also obvious that the absolute inaccuracy A^cm of the Moliere's 
theory in determining the screening angle Xa with the rise of Z increases to 
20% and the corresponding relative inaccuracy 6ccm is about 34 percent. 
Thus, in the case of scattering on targets with large Z the such corrections 
as AccM (Xa) and 6ccm become significant and should be taken into account 
in determining the lifetimes of relativistic elementary atoms in experiments 
with nuclear targets. 

4 Conclusion 

1. Within the framework of unitary Glauber approximation for particle- 
atom scattering, we develop the general formalism of the Moliere-Fano 
multiple scattering theory. 

2. We have estimated the relative unitary corrections to the parameters 
of the M-F theory, resulting from a reconstruction of it's unitarity in 
second-order optical-model of the Glauber theory and found that they 
are of order Za^. 

3. We have also received the rigorous relations for the exact and the Born 
values of the quantities Qeiijf)^ in [(7(77)], and xLi which do not depend 
on the shape of the electron density distribution in the atom and are 
valid for any atomic model. 

4. We have calculated the Coulomb corrections A^c = A^c [in (Xa)] — 
Ace [Infi-lXc)] and relative corrections 5cc = Scc{xL) = ^cc [giXc)] = 
Sec i^MiXct)] with nuclear charge ranged from Z = 2 to Z = 82 and 
showed that for Z ~ 80 these corrections comprise the order of 30% 
and 40%, correspondingly. 

5. Additionally, we evaluated the absolute and relative inaccuracies of the 
Moliere's theory in determining the screening angle and found that for 
heavy atoms of the target material {Z ~ 80) they are about 20% and 
30%, respectively. 
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